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ABSTRACT 



Context. The molecular species hydrogen peroxide, HOOH, is likely to be a key ingredient in the oxygen and water chemistry in the 
interstellar medium. 

Aims. Our aim with this investigation is to determine how abundant HOOH is in the cloud core p Oph A. 

Methods. By observing several transitions of HOOH in the (sub)millimeter regime we seek to identify the molecule and also to 
determine the excitation conditions through a multilevel excitation analysis. 

Results. We have detected three spectral lines toward the SMI position of p Oph A at velocity-corrected frequencies that coincide 
very closely with those measured from laboratory spectroscopy of HOOH. A fourth line was detected at the 4cr level. We also found 
through mapping observations that the HOOH emission extends (about 0.05 pc) over the densest part of the p Oph A cloud core. We 
derive an abundance of HOOH relative to that of H2 in the SMI core of about 1 x 10"'". 
Conclusions. To our knowledge, this is the first reported detection of HOOH in the interstellar medium. 
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1. Introduction 

Hydrogen peroxide, HOOH, is believed to play an important 
role in the Earth's atmospheric ozone and water chemistry. It 
is a key constituent in the gas- and liquid-phase radical chem- 
istry and has an oxidizing potential in the liquid phase. Gas- 
phase HOOH has been se en in the Martian atmosphere by 
groun d-based observations (IClancv et al.l l2004t lEncrenaz et alj 
I2OO4I) . Howe ver, recent Mars obs ervations with the Herschel 
Observatory (Hartogh et"aD 1201 Ol) failed to detect HOOH at 
levels below those previously seen. The non-detection was at- 
tributed to seasonal variations. 

Interestingly, HOOH is among the simplest molecules that 
show internal rotation. The internal rotation, or torsion, mani- 
fests itself as a rotation of the two 0-H bonds about the O-O 
bond. This hindered internal rotation can be desciibed with a 
torsion potential in which the two minima (the most stable con- 
figurations) do not coincide with the cis or trans alignment of 
the two O-H bond^. The twofold barrier gives rise to a quartet 
of sublevels fo r each torsional sta te. These sublevels are denoted 
T = 1,2,3,4 jHunt et al.iri965h . Moreover, HOOH is a light, 
slightly as ymmetric prola te rotor (belonging to the C)^^ point 
group, see Housed (11984 for a discussion), with only c-type 
transitions and with a dipole moment of 1.6 D (iCohen & PickettI 



[T98TI iPerrinet al.l[T996l) . Because of this, no pure rotational 
transitions occur, and only transitions corresponding to a com- 
bined ro-torsional motion change can take place. The mm and 
submm spectrum of HOOH has been studied rather extensively 
in the laboratoiy (iHelminger et al . 1981; Petkie et al.l l 19951) and 
is available in the JPL database (iPickett et aL.,1998) . 

To our knowledge, HOOH has not so far been detected in the 
interstellar medium. Also, very few abundance limits have been 
reported. Blake et al. (1987) derived an upper limit for HOOH 



of 4.5 X 10 with respe ct to H2 from their Orion spectral scan 
. iBoudin et al.l d 19981) reported an upper limit of 5.2 % of 



data. 



* Based on observations with the Atacama Pathfinder Experiment 
(APEX) telescope. APEX is a collaboration between the Max-Planck- 
Institut fiir Radioastronomie, the European Southern Observatoiy, and 
the Onsala Space Observatory. 

' When the two O-H bonds point in the same direction, this is refen'ed 
to as the cis position, while the 180 degree opposite case is called the 
trans position. For HOOH, the trans potential barrierheight is 557 K, 
while the cis barrier height is almost 4000 K dPelz et al.il993h . 



solid HOOH relative to H2O ice toward NGC 7538 IRS9. 

As in the Earth's atmosphere, HOOH is expected to be 
closely connected to the water and molecular oxygen chemistry 
also for those physical conditions prevailing in molecular clouds. 
In cuiTent pure gas-phase models, HOOH is formed by reaction 
of H2 with HO2 or via the reaction involving two OH radicals. 
Howeve r, these r eactions proceed very slowljB Alternatively, 
Tielens & HagenI (11982) suggested that HOOH could be formed 
on grain surfaces by the successive additions of H atoms to 
molecular oxygen. If this is the case, HOOH could be closely 
related to the amount of O2 on grains. 

The 119 GHz line of O2 was detected toward the cloud 
p Oph A with an a bundance of 5 x 10 relative to H2 by 
iLarsson et al. I (l2007l) using the Odin satellite. The p Oph A 
molecular cloud, at a distance of about 120 pc, has been the 
subject of se veral studies. Continuum observations (Andre et al 
1993t iMotteet al.lll998h and C"^0 observations (Lise au et al 
201(|) revealed several cores. Very recently, Bergman and 
coworkers found a very high degree of deuteration toward the 
SMI core in p Oph A from observations of deuterated H2CO 
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Table 1. Observed HOOH lines 



Frequency 




Transition 






E„ 


A,i 




(MHz) 




- J"„ „ 


t' - 


r" 


(K) 


(s- 


) 




219166.86 


3o,3 


-2i,i 


4- 


-2 


31.2 


8.58 X 


10- 


-5 


251914.68 


6l,5 


~ 5o,5 


2- 


-4 


65.5 


2.46 X 


10- 


-4 


268961.17 


4o,4 


-3l,2 


4- 


-2 


41.1 


1.84 X 


10- 


-4 


318222.52 


5o,5 


-4l,3 


4- 


-2 


53.4 


3.31 X 


10- 


-4 


318712.10 


5 1,4 


- 6o,6 


3- 


- 1 


67.0 


4.12 X 


10- 


-4 


670595.82 


li.o 


-Oo.o 


3- 


- 1 


32.2 


5.79 X 


10- 


-3 



(iBergm an et alj 1201 lb . These observations suggest that grain 
surface reactions were at work to produce the very high deu- 
terium levels observed in the gas-phase material. 

In this Letter we continue our study of the p Oph A cloud by 
reporting on observations of several HOOH transitions. We have 
detected four HOOH lines and two of these lines were mapped 
over the central part of this cloud. In Sect. ([2} we describe our 
observations and present our results. There, we also describe in 
more detail the energy level structure and symmetry of HOOH, 
which is important for discussing the detected lines. Then, in 
Sect. (O, we discuss the implications of our results. 



2. Observations and results 

We have used the APEX 12 m telescope located at abou t 
5100m altitude in the Chilean Andes (Gusten et al., |2006|) 
to observe HOOH. For the lower frequency Unes we used 
the Swe dish heterodyne fac ility instruments APEX-1 and 
APEX-2 jVassilevet al.ll2008l) . The 7-pixel longer wavelength 
(450 jum) module of the MPIfR-built CHAMP+ receiver array 
( Kasemann et alJl2006) was used for the observations of a high- 
frequency line. The 7 pixels, spaced by 18", are arranged in a 
hexagon around a central pixel. The observations took place on 
several occasions during 2010; April 2-11, July 7, August 4-8, 
and September 10-13. 

The targeted HOOH lines are listed in Table [T] The fre- 
quency, transition quantum number designation, energy of upper 
level, and Einstein A-coefficient are listed. All values have been 
compiled from the JPL catalogue (Pickett et al. 1998). The fre- 
quency uncertainty for the listed lines is 0. 1 MHz or better and 
this corresponds to 0.14kms ' at 219 GHz. We also present in 
Fig. [T] the energy diagrams for all levels below 100 K. Owing 
to the symmetry of HOOH, four different radiatively decoupled 
ladders occur: Ai„3, A2m4, Bim3, and B2«4- The subscript in- 
dicates the pair of torsional quantum numbers t involved. The 
c-type elec tric dipole tr ansitions must also obey r = 1 <-> 3 or 
T = 2 4 (iHunt et al.ll 965) and are drawn as downward arrows 
connecting the upper and lower levels in Fig. [T] The A-species 
have a nucl ear spin weigh t of 1 and the B-species have a spin 
weight of 3 (lHougenll984l) . This is, of course, due to the nuclear 
spin directions of the two H atoms (in the same way as the ortho 
and para symmetries occur for H2O or H2CO). 

The ground-state symmetry species is that of A 1^3 with the 
Bi„3 state only about 2.5 K higher energy. This means that the 
energy difference between the species with different nuclear spin 
weights is much smaller than the corresponding difference for 
H2O or H2CO. The torsional r = 2,4 states are about 16 K 
above the t = 1, 3 states. This difference stems from the tunnel- 
ing through the trans barrier and is comparabl e to the gas kinetic 
temperatures of 20 - 30 K found in p O ph A dLoren et al.lll990t 
iLiseau et alJl2003HBergman et al.ll20Tlh . 



HOOH energy diagram, spin weight=1 




1 2 ff^ 1 2 

1 3 1 r 4 2 4 
HOOH energy diagram, spin weight=3 




1 2 jr^ 1 2 

1 3 1 T 4 2 4 



Fig. 1. HOOH energy level diagrams. The energy is given in K 
on the vertical axis, and at the bottom the quantum numbers Kg 
and T are shown. To the left of each level the rotational quantum 
numbers Jk„,k, are listed. The upper diagram shows the levels 
with a nuclear spin weight of 1 (A 1^3 and A2M4), while those 
levels with a spin weight of 3 (Bim3 and B2«4) are shown in 
the lower diagram. The allowed c-type transitions within each 
of the four ladders are indicated as arrows. The red arrows rep- 
resent detected lines, while the red-dashed arrows indicate non- 
detections. 



In Fig. |2] we show the HOOH spectra toward the core of 
p Oph A. The upper five spectra are toward the SMI core, 
Q'(J2000) = le*" 26" 27.2' and 5(J2000) = -24° 24' 04". The 
670 GHz CHAMP^ spectrum is an average of all pixels and 
is centered 30" north of the SMI position (usually denoted 
SMIN). In Fig. [3] spectra of the 219166 MHz line for the ob- 
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5„,-4,, 4-2 



219.166 GHz 




251.914 GHz 




318.222 GHz 



5, .-6„. 3-1 



318.711 GHz 



4„,-3,^4-2 




2B8.960 GHz 



l,„-0„„ 3-1 



670.592 GHz 



^LSR [km s ] 

Fig. 2. HOOH spectra toward p Oph A. The transition is in- 
dicated in each spectrum. The intensity scale is in K and 
the velocity (vlsr) scale is in kms"'. The velocity resolution 
is 0.25 km s"'. 

Table 2. Observed HOOH line velocities, widths, and intensities 



Freq. 


Beam size 


Vlsr 


FWHM 




(MHz) 


(arcsec) 


(kms-') 


(kms-') 


(Kkms-') 


219167 


28 


3.8 


0.84 


0.167(0.018) 


251915 


25 


3.7 


0.75 


0.165(0.018) 


268961 


23 


3.7 


(1.2) 


0.040(0.011) 


318223 


20 


3.8 


0.78 


0.106(0.013) 


318712 


20 






< 0.044 


670596 


9 






< 0.12* 
<0.05r 



^average value from three pixels closest to the SMI position 
''average value from all pixels, central pixel on SMIN 



served map positions are shown. In the southern part the line is 
centered at 3.7 km s"', while further north, at offset (0", +60"), 
the line is at 3.2 km s '. This NS velocity gradient is almost 
identical to the one s een for H2CO and its deuterated variants 
(^Bergman et al."201 1). The narrow peaks to the NW are adjacent 
to where the sulphur species peak as noted by the same authors. 
In addition to the map in Fig. [3] we also mapped the 251 GHz 
line, albeit with a poorer S/N ratio and will not discuss it further 
here. 

The beam sizes, velocity-integrated main-beam brightness 
temperatures (J T^^ dv) , fitted LSR velocities (vlsr), and line 
widths (FWHM) are tabulated in Table |2l The r^b-scale was 
established assuming main-beam efficiencies of 0.75, 0.73, and 
0.4, for APEX-1, 2, and CHAMP+, respectively. The upper lim- 
its are 3cr and for the CHAMP^ line two intensities are listed in 
Table |2j one for the three pixels closest to the SMI position and 
the other by averaging data from all pixels. The tabulated errors 
and upper limits depend only on the channel noise. 
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Fig. 3. Map spectra of the 219166 MHz HOOH line toward 
p Oph A. The map offset (0", 0") corresponds to the SMI posi- 
tion. The intensity scale, in K, and the velocity (vlsr) scale, 
in km s are indicated in the upper right spectrum. The velocity 
resolution is 0.25 km s"'. 



Using the integrated line intensities in Table |2] (corrected by 
a beam-filling factor corresponding to a source size of 24", see 
Bergman et al l (1201 1^1. we perfo rmed a rotation diagram anal- 
ysis (Goldsmith & Langei|ll999|). From this we can determine 
the rotation temperature, T^ot, as well as the HOOH column den- 
sity, A^(HOOH). The resulting rotation diagram is displayed in 
Fig. |4] The fit is based on the detected t = 2 4 lines. 
We derive T^ot = 22 + 3 K and a total HOOH column density 
of (8 + 3) X 10'^ cm"^ where the errors depend on the uncer- 
tainty of the integrated intensities and a calibration uncertainty 
of 10 - 15%. This rotation temperature is very close to the esti- 
mates of the kinetic temperatures from H2CO at the SMI posi- 
tion. Obviously, the non-detection of the t = 1 3 lines (open 
squares in Fig. HJi is not consistent with the fit. Especially the 
670 GHz line should have been detected given these values of 
r,ot and A^(HOOH). Hence, we conclude that the r = 1, 3 states 
are not populated according to a simple LTE model. The non- 
detection of the 670 GHz line could be a result of subthermal 
excitation. The spontaneous rate coefficient is large (TablelTJ and 
for a typical value of a collision coefficient (~ 10"'° cm^ s"'), the 
critical density of the 670 GHz transition is about two orders of 
magnitude higher than the H2 density determined from H2CO 
and CH3OH observations in the same source (Bergman et^ 
2011). From the energy diagrams (Fig. [T) it is clear that there 
is a lack of radiative de-excitation routes out of the Kg = Q levels 
in the t = 1, 3 states (as opposed to the t = 2, 4 states). The pop- 
ulation of the T = 1 , 3 states could therefore be confined to the 
K„ - levels if the collisional excitation is inefficient. However, 
a full statistical equilibrium analysis is needed for understand- 
ing the details of the excitation. This also requires some basic 
knowledge of the collision coefficients. For now we assume that 
the T = 1,3 population is negligible compared to the r = 2,4 
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; Total N(HOOH) = (7.9 ±2.8) 10'^ cm"^ E 
N(HOOH, T=2-4) = (2.5 ±0.9) lO'^cm"^ E 

20 30 40 50 60 70 

Upper energy (K) 

Fig. 4. HOOH rotation diagram. The detected lines (only t - 
2 4 lines) are shown as open triangles with error bars. The 
size of the error bars corresponds to the total uncertainty (noise 
and calibration). The 3cr upper limits (not included in the fit) are 
shown as open squares and downward arrows and originate from 
transitions between t - 1,3 states. The fitted rotation tempera- 
ture is noted, as is the total column density. The column density 
for unpopulated t = 1 , 3 states is also given. 



population. In this case, we instead find a total molecular col- 
umn density of (3 + 1) x 10'^ cm"^. 

3. Discussion 

Given the good agreement of the velocities of our detected 
HOOH lines (together with the w 0. 1 km s ' accuracy of the lab- 
oratory frequencies) with those from other species we are very 
confident that the lines belong to HOOH. From our mapping of 
the 219 GHz line it is also evident that the NS velocity gradient 
seen for HOOH reflects that of other species. Moreover, the de- 
rived rotation temperature of 22 + 3 K is what would be expected 
for the p Oph A cloud core. For this fairly low excitation temper- 
ature one would not expect many lines from other speci es to be 
present and, using the JPL and Cologne databases ( Picke ttet al.l 
\199^ : Muller et al. 2005), we found no lines from other species 
that could possibly interfere with the identification. Of course, 
the narrow lines (with FWHM typically < Ikms ') seen to- 
ward the p Oph A cloud core also make line confusion much 
less likely. 

From the H2CO and CH3OH analysis of the SMI core 
iBergman et al.l (l201lb determined an H2 column density of 3 x 
lO^^cm"'^. Assuming that the HOOH level populations mainly 
reside in the r = 2, 4 states, we arrive at an HOOH abundance of 
about 1 X 10""'. This is well below the l imit of 4.5 x 10"'° found 
toward Orion KL bv lBlake et all (Il987h . 

According to current gas-phase schemes (e.g., the OSU 
chemical reaction database), formation of HOOH in the gas 
phase is not efficient. Only two very slow reactions are proposed 
for its formation, reaction of H2 with HO2, or reaction of two 
OH radicals. 

On grains, HOOH is formed through successive hydrogen 
additions to O2. This was first proposed by Tielens & HagerJ 
([1982), 

based on theoretical arguments. Recent laboratory exper- 
iments have been made to investigate this route, up to formation 
of water molecules: 

O2 -H H -H H ^ HOOH (1) 
HOOH -H H ^ H2O + OH. (2) 

iMivauchi et al.l (l2008l) have investigat ed the reaction of H 
atoms with solid O2 at 10 K. Subsequentlv. lloppolo et alJ ( l2008l) 



have investigated the same reaction in the temperature range 12- 
28 K. Both studies showed that the conversion of O2 stops at 
some point before exhaustion of O2 beca use of shielding of O2 
in the deepest layers. The experiment of lloppolo et al.l (l2008h 
shows that the shielding decreases with increasing temperature, 
pointing to the fact that the O2 ice may become more porous 
when close to its sublimation temperature (30 K). 

Because this shielding may not be very relevant in space, 
lOba etalJ(l2009l) investigated the formation of HOOH and H2O 
when codeposing O2 and H in the temperature range 10-40 K. 
The H2O/HOOH ratio in the formed ices is observed to depend 
strongly on the temperature, and on the O2/H flux. The mea- 
sured H2O/HOOH ratio is lower than 5 in all experiments {T = 
10, 20 K and 02/H-flux between 3.8 x lO""* and 1.9 x 10"^), al- 
though higher values may be obtained in the case of lower 02/H- 
flux. The present detection of HOOH may open the possibility to 
quantify the importance of reactions (1) and (2) in the formation 
of water. 

lopp olo et alJ (l2008l ) have modeled the formation of water in 
typical dense clouds (their Figure 4). The laboratory results lead 
to a revision of the energy barriers involved in the models, and 
their new model (only accounting for the three main routes of 
water formation on the grains) predicts a fractional abundance 
for HOOH of a few 10"'^ with respect to H nuclei. This is more 
than three orders of magnitude lower than our detection. 

Further understanding will require detailed chemical model- 
ing of grain chemistry. This will be the scope of a forthcoming 
paper. The observation of water in the p Oph A region with the 
Herschel Observatory as well as the confirmation of the O2 de- 
tection would also be very valuable in setting constraints on the 
models. 
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